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Abstract. We discuss the possibility to constrain the relation between redshift and temperature of the cosmic 
microwave background (CMB) using multifrequency Sunyaev-Zeldovich (SZ) observations. We have simulated a 
catalog of clusters of galaxies detected through their SZ signature assuming the sensitivities that will be achieved 
by the Planck satellite at 100, 143 and 353 GHz, taking into account the instrumental noise and the contamination 
from the Cosmic Infrared Background and from unresolved radiosources. We have parametrized the cosmological 



temperature-redshift law as T oc (1 + z 



Jl-a) 



Using two sets of SZ flux density ratios (100/143 GHz, which is 
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most sensitive to the parametrization of the T — z law, and 143/353 GHz, which is most sensitive to the peculiar 
velocities of the clusters) we show that it is possible to recover the T — z law assuming that the temperatures and 
redshifts of the clusters are known. From a simulated catalog of ~ 1200 clusters, the parameter a can be recovered 
to an accuracy of 10 -2 . Sensitive SZ observations thus appear as a potentially useful tool to test the standard 
law. Most cosmological models predict a linear variation of the CMB temperature with redshift. The discovery of 
an alternative law would have profound implications on the cosmological model, implying creation of energy in a 
manner that would still maintain the black-body shape of the CMB spectrum at redshift zero. 



Key words, cosmology: 
microwave background 



large-scale structure of Universe, galaxies: clusters: general, cosmology: cosmological 



1. Introduction 

One of the strongest predictions of the standard cos- 
mological model is the linear variation of the Cosmic 
Microwave Background (CMB) temperature with redshift: 
Tcmb{z) — Tq(1 + z) where To is the temperature of 
the CMB at redshift zero. T> has been determined with 
a high accuracy from data taken by the Far Infrared 
Absolute Spectrometer (FIR AS) onboard the COB E satel- 
lite: T = 2.725 ± 0.002 K l|Mather et al.lll999|) . Those 
measurements do not imply that the CMB spectrum is 
Planckian at all redshifts, nor give information about how 
the temperature varies with redshift. Alternative cosmo- 
logical models imply a different redshift dependence of the 
CMB temperature, for instance: 



r (i + z) 



(i- 



with < a < 1 



(1) 



IlLima et al . 2000). Equation Q requires energy injection 
into th e CMB, which cou ld come from decay of vacuum 
energy l|Freese et al.lll987|) . In such a case, the Planckian 
form of the radiation spectrum is preserved if the creation 
of photons is adiabatic, i n the sense t hat the entropy per 
photon remains constant jLimalll99dh 
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An important test of such models would be to deter- 
mine the temperature of the CMB at different redshifts 
through observations. 

Shortly after the discovery of the CMB by 
IPenzias fc Wilsonl l)l965[) it was realized that the relic ra- 
diation had already been indirectly observed through ob- 
servations of interstellar abso rption lines with very low 
excitation temperatures (see iThaddeusl Il972l for a re- 
view of early observations). More recently, a temperature 
T = 2.729^0^31^ has been derived from a survey to de- 
termine the a mount of CN rotat i onal excitation in inter- 
stellar clouds fe.oth et alJll99l . iBahcall fc Woll lll9fifift 
pointed out that the temperature of the CMB could be 
probed at much higher redshifts by studying the absorp- 
tion spectra of quasars and searching for the species with 
energy states whose populations are sensitive to the am- 
bient radiation temperature. This method, however, gives 
only an upper limit on the CMB temperature unless the 
local contribution to the excitation can be estimated, 
since other processes can contribute to the excitation of 
the atoms /molecule s llGe et al.ll997llSrianand et all2 000. 
iMolaro et all [20021. So far, the standard law is consis- 
tent with observations, but so are alternative laws with 
a = 0.003 ± 0.13 or b = 0.99 ± 0.22 for a s traight lin e 
T(z) = T (l + bz) l|Losecco et al J200li see also lPuvl200^) . 
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Central frequency (GHz) 
Bandwidth (GHz) 
Angular resolution (arcmin) 
Effective beam area (arcmin 2 ) 



fOO 143 217 353 

33 47 72 116 

9.2 7.1 5.0 5.0 

95.9 57.1 28.3 28.3 



27.0 



Random noise (mjy/beam) 14.0 10.2 14.3 

Cosmic Infrared Background (mjy/beam)10.5 15.6 23.7 81.9 |jf 

Unresolved radiosources (mjy/beam) 26.0 15.5 8.6 - ^ 

Total uncertainty (mjy/beam) 31.4 24.3 29.0 86.3 I 

f 

Table 1. Four frequency channels of the Planck high- < 
frequency instrument relevant to this paper. The effective 
beam area has been calculated assuming a Gaussian beam. 
The levels of contamination due to the Cosmic Infrared 
Bac kground and to un resolved radiosources are taken from 
lAghanim etaH l)2004j) . 



Another way to probe the temperature of the 
CMB is provided by the Sunyaev-Zeldovich (SZ) effec t 
i Zeldovich fc Sunvaevlll969l ISunvaev fc Zeldovichlll97l . 
The CMB spectrum is distorted due to interaction with 
intervening hot electrons. Massive clusters of galaxies con- 
tain significant amounts of hot gas that can produce a 
measurable SZ signal (see Fig. Both the decrement 
at centimeter wavelengths and the increment in the mil- 
limeter range have been detected and large dedicated SZ 
surveys are planned (s ee lRephaelilll995l lBirkinshawlll999t 
ICarlstrom et al.l 120021 for reviews). Using the standard 
scaling v = isq(1 + z) and the temperature-redshift relation 
given in eq. ljT|l. the dimensionless frequency x = hv/kT 
depends on z, and we have 



kT ~ kT Q [ + ' ' 



(2) 



where h is the Planck constant and k is the Boltzmann 
constant. The redshift independence of the SZ effect there- 
fore occurs only when the standard temperature law is 
assumed. 

Two methods can provide information on the CMB 
temperature at the redshift of t he intervening cluster 
l|Fabbri et alJll978llRephaelilll980|) . The first one involves 
a precise measurement of the null point of the SZ signal, 
located around 218 GHz. The precise frequency depends 
both on the temperature of the hot gas and its redshift 
(see Sect. It also depends on the peculiar velocity of 
the cluster. The second method involves measurement of 
the SZ signal at several frequencies; the temperature of the 
CMB can be derived from the ratio of the observed SZ flux 
densities in different bands, which for some frequencies de- 
pends weakly on the properties of the intracluster gas (see 
Sect. [P. 

In this paper, we investigate the possibility of apply- 
ing these methods to upcoming SZ observations such as 
those planned in th e Planck mission (see Table 1 , and 
lLamarre et all l200.i iTauberl Eiol iMennella et "all 120041 
about this mission). 
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Fig. 1. Spectral signature of both the thermal SZ effect 
(solid line) and the kinetic effect (dashed line) for the 
standard cosmological temperature law (a = 0). The dash- 
dotted and the dotted lines show the corresponding curves 
for a = 0.1 and z = 0.5. The Planck frequencies used in 
this study are indicated. The following cluster parameters 
were used: kT e = 5 keV, v r = +1000 km s _1 . 

2. The SZ effect 

The thermal SZ effect can be described in terms of the 
Compton y parameter: 



y = a T 



kT K 



m e c* 



dl 



(3) 



where <jt — 6.65 x 10~ 25 cm 2 is the Thomson cross section, 
n e is the number density of electrons, T e is the tempera- 
ture of the intracluster gas, m e is the electron mass, and 
c is the speed of light. The y parameter is a measure of 
the integrated pressure of the intracluster gas along the 
line-of-sight. 

Due to inverse Compton scattering off the hot elec- 
trons, the spectrum of the CMB is distorted. In the non- 
relativistic limit (if the velocity of the hot electrons is neg- 
ligible with respect to c), the change in the CMB flux 
density is 



S v (x) = S^ a {x)Q{x)Y 



(4) 



where S™ B = 2{kT) 3 x 3 /{{hc) 2 (e x - 1)) is the unper- 
turbed CMB flux density. The spectral distortion of the 
incident CMB spectrum is described by the function 



Q(x) 



1 \tanh(a;/2) 



-4 



(5) 



Y is the integrated Compton parameter over the projected 
area of the cluster on the sky: 



y = d: 



yd A . 



(6) 



where Da is the angular diameter distance. 

For high electron temperatures the relativist ic effect 
has to be taken into account (e.g., IWrightjri979t iFabbril 
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Fig. 2. Variation of the cross-over frequency as a func- 
tion of the redshift of the intervening cluster for different 
cluster temperatures and different T — z relations (the pa- 
rameter a varies from to 0.1 by 0.02 from top to bottom). 
Here the peculiar velocity of the clusters has been taken 
as equal to 0. 
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Fig. 3. Variation of the cross-over frequency as a func- 
tion of the peculiar velocity of the intervening cluster for 
different cluster temperatures and different T — z relations 
(the parameter a varies from to 0.1 by 0.02 from top to 
bottom). 



ll98lllRenhaeli fc Yankovit chill 9971) . Useful numerical ta- 
bles and analytic fit ting formulae have be en provided by 
several authors fsee lltoh fc Nozawall2004l and references 
therein). 

The kinetic SZ effect is due to the motion of the cluster 
with respect to the CMB frame: 



AT _ v r 
T c ' 

with the convention that v r 



(7) 



> for a receding cluster. 



3. Measurement of the cross-over frequency 

The location of the null point of the SZ signal v cr can give 
information on the temperature of the CMB at the redshift 
of the intervening cluster if it doesn't obey the standard 
law. From cq. (GJ we have 



r*^(l + *)- 



a . Sazanov 



& Sunyaev ( 1998a ul998bl ) provide a simple approximation 
for the non-dimensional cross-over frequency x c 



3.830 1 + 0.30 



v r m e c 
c kT e 



1.1- 



kT e 
m e c 2 



1.5^ ) (8) 



where v r is the velocity of the cluster along the line of 
sight and T e is the temperature of the intraclus ter gas. 
For v r = one recovers t he approximation given bv lFabbril 
(Il98lh . lltoh et alJ ( 1998) provide a second-order fit to the 
numerical result, x cr = 3.83 (l + 1.16746> e - 0.85336^) , 
valid for < kT e < 50 keV and more accurate at high T e 
than eq. and where 9 e = kT e /m e c 2 . In the following 
we use that approximation in the figures where v r = 0. 



Figure El shows the variation of the cross-over fre- 
quency as a function of the redshift of the intervening clus- 
ter for different cluster temperatures and different T — z 
relations. Here we first consider v r = 0. For a — 0, one re- 
covers the known result that v cr is constant with redshift 
and increases with increasing T e . For different values of a, 
the curves are split and v ct decreases with increasing a. 
There is thus a parameter degeneracy between T e and a. 

FigureElillustrates the dependence on the peculiar ve- 
locity of the cluster, v r , introducing a further degeneracy. 

Those parameter degeneracies and the large band- 
widths of current bolometers make it difficult to apply this 
method in practice to constrain the CMB temperature. 

4. Measurements at different frequencies 

A more promising method consists of using ratios of ob- 
served flux densities in several frequency bands. It has 
recently been applied sucessfully to determine the CMB 
temperature at the redshift of the Com a cluster (at z = 
0231 ) and of Abell 2163 at z = 0.203 bv lBattistelli et alJ 
(2002) who used observations at 32, 143 and 272 GHz. 

We have calculated the flux density of the thermal 
SZ effect using the analytic fitting formulae given by 
lltoh et aD (^998) valid in the relativistic case. We have 
added the contribution of the kinetic effect. The flux den- 
sity ratio depends weakly on the cluster temperature for 
frequencies between 32 GHz and ~ 170 GHz. If the second 
measurement (1/2) is done at a higher frequency, the flux 
density ratio decreases with increasing T e , more strongly 
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Fig. 4. Variation of the flux density ratio as a function 
of the cluster temperature for clusters at z = 0.5 and 
different z — T relations. The parameter a varies from 
to 0.1 by 0.02 from bottom to top. Here the peculiar ve- 
locities have been set to zero. The dashed lines show the 
100/143 GHz flux density ratio, multiplied by —0.6 for 
clarity. The dash-dotted lines show the 143/353 GHz ra- 
tio. The dotted lines show the 100/353 GHz ratio. Clearly, 
the 100/143 GHz ratio is most sensitive to a and varies 
least with the cluster temperatures. 



as i/2 increases. For different values of the a parameter, the 
curves are split. The magnitude of the split increases with 
redshift. Figure 0] shows the variation of the flux density 
ratios with the cluster temperature, Fig. the variation 
with the peculiar velocity of the cluster. Observations with 
new X-ray satellites (XMM, Chandra) make it possible 
to measure T e with good accuracy. Note, however, that 
the flux density ratios depend only weakly on the cluster 
temperature. To constrain the parameter a, it would be 
necessary to constrain the peculiar velocity of a cluster 
independently. Alternatively, it is possible to constrain a 
by observing a large number of clusters with known gas 
temperatures, as will be shown in the next section. 



5. Simulating a large SZ survey 

We seek to simulate the yield of a Planck-like SZ survey to 
test how effectively our parametrization of the T — z rela- 
tion can be constrained. One way to do this is to simulate 
sky maps based on iV-body simulations with SZ signals as 
well as random and systematic noise in the relevant fre- 
quency bands and then use a suitable cluster extr action 
algorithm (see, for example. iGeisbiisch et al]l2005|) . Here 
we take the simpler approach of directly simulating a clus- 
ter catalog from the redshift and mass distribution of dark 
matter halos and computing the SZ signals assuming that 
the clusters are unresolved. This method is justified for 
this study as it gives approximately the same number of 
recovered clusters as more general methods. 
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Fig. 5. Same as Fig.0]but as a function of the peculiar ve- 
locity of the cluster. The 100/143 GHz flux density ratios 
have been multiplied by 0.5 for clarity. Here kT e = 5 keV. 
The 143/353 GHz flux density ratio is most sensitive to 
the peculiar velocities. 



5.1. The model 

We have used a fiducial ACDM cosmology with Q = 
0.3, A = 0.7, H = lOOh km s^Mpc" 1 , h = 0.7, 
i ts = 0.9, consisten t with the first-yea r WMAP results 
jSnergel et alJl2003h . We have used the ISheth fc Tormenl 
(1999) cluster mass function, which provides an excel- 
lent match to the results from large numerical simula- 
tions, and has the advantage of having a theoretical jus- 
tification in terms of the ellipsoidal collapse of clusters 
llSheth et all 1200 ll). The mass function can be expressed 
as llMo fc WhiteJl2002lk 



it M dM 



(9) 



where po is the comoving background density, v 1 — ^/av, 
a = 0.707, A ~ 0.322, q = 0.3, v = v(M,z) = 
5 c /(D(z)<t(M)), S c is the linear overdensity of a pertur- 
bation at the time of collapse and virialization; S c varies 
weakly with the cosmological parameters and has a value 
of 8 C ~ 1.69. a(M,z) is the variance of the fluctuation 
spectrum on a mass scale M = |7ri? 3 p - 

The number density of clusters above a certain mass 
per unit steradian per redshift interval is 



d 2 N(z) d 2 V c 



dildz dQdz 



-foe 



n(M, z)dM 



Jlim(z) 

where the comoving volume is 



d 2 v c _ D 2 A ( z )(i + z y 



dQdz 



H(z) 



(10) 



(11) 



where Da is the angular diameter distance and H(z) is 
the Hubble parameter. 
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Fig. 6. Left: Limiting mass as a function of redshift at three frequencies for observations at 3 a level using Planck, 
taking into account both the random instrumental noise and the levels of contamination due to the Cosmic Infrared 
Background and to unresolved radiosources given in Table 1. The thin lines correspond to the relativistic calculation 
and the thick lines to the non-relativistic one. 

Right: Number density of clusters per redshift interval and per square degree on the sky as a function of redshift, above the 
limiting mass. 



A mass-temperature relation is needed in order to as- 
sign a gas temperature to a cluster. If the gas is in hydro- 
static equilibrium with the gravitational potential of the 
cluster, then: 



kT e 



1 



din p gas ( r )/d\nr\ rvil 



fimp- 



GM V 



(12) 



where /3 gas is the gas density, r v j r is the virial radius, and 
[i = 4/ (5 A + 3) is the mean molecular weight, X is the 
hydrogen mass fraction a nd mp is the mass of the proton. 

We have adopted the iNavarro et al' ] lll997h profile for 
the gas: 



din 



Pgas 



d\nr 
with 



l + 3c 
1 + c 



M vir 



10 14 /l- X Mf; 



-1/5 



(13) 



(14) 



d In p g 
On? 



Taking a constant c = 5 yields an error in 
eq. (|T3|> of less than 5% for the range of masses considered. 

Introducing A c , the average density within the cluster 
relative to the critical density of the background at the 
redshift of the cluster, one obtains: 



kcV 



= 2| 

( 6.8 



7.75 



\^5X+3 

x(l + z) 



d 111 /9 gas (r) I d In ^|r- vir 

( Mv, r 1 



2/3 



(15) 



TUT) 



1/3 /A^ 1/3 
V 178/ 



We have cal culated A r from Ay , ,. = p i pbg, taken from the 
fits given bv lKitavama fc Sutol (H.996): 



18t^(1 + 0.4093?< r yu ^) , 



(16) 



1, and f2 v i r is the density parameter 



Oo(l + Zvirr 
fi (l + Zvir) 3 + A 



for a flat ACDM 



where w V i r = l/il v ; ; 
at virialisation; fi v 
model. 

We have added a spread of 30% to the mass- 
temperature relation. 

As for the peculiar velocities of the clusters, we have 
drawn them from a Gau ssian distribution with zero mean. 
ISheth & Diaferiol ( 200l]) provide a fit for the dispersion as 
a function of the mass of the cluster: 



<7Wo(M) = 



Cfit 



1 + {R/R 



lit J 



(17) 



with R = (3M/47rpo) 1/3 , cr fit = 414.7 km s _ \ Rm = 
34.67/1" 1 Mpc, and r\ = 0.85 for a ACDM model with 
the same parameters as ours. cr^ alo is the dispersion of 
the Maxwellian distribution (three-dimensional peculiar 
velocity distribution). For the peculiar velocity along the 
line-of-sight, we thus took <7p ec = ^ alo /3. 

5.2. The limiting mass 

The integration of eq. (|10|) is to be done from a limit- 
ing mass which is redshift-dependcnt. The majority of the 
clusters will be unresolved at the resolution of the Planck 
instruments l(Kav et alJl200ll) . 
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Fig. 7. Flux density ratio versus redshift for simulations with different values of the parameter a. For clarity, only 
instrumental noise has been included. The triangles correspond to a — 0, the plusses to a — 0.1, the crosses to a = 0.2. 
Left: In the absence of peculiar motions (v r = 0). The distributions are clearly different. Right: Including the peculiar 
motions. The spread in the flux density ratio is much larger than when v r = 0, and the different a distributions are 
mixed at low redshift s. 



The SZ flux can be expressed as 

Sv = 2.29 x 10 4 -£Lq(.t) x 1.70 x 10~ 2 /i ( 

/ n-2 V ' (18) 

x (m{m^) (i^-Jy, 

where /icm is the gas mass fraction that we will take con- 
stant and equal to 0.1 and kT e is taken from eq. (|15fl 
l|Fan fc Chiuehll200lj) . 

We have calculated the limiting mass for Planck ob- 
servations at a detection level of 3 a, where a is the total 
uncertainty given in Table 1. The error estimates on the 
flux densities are discussed in Sect. 15.41 

5.3. The SZ catalog 

We have distributed clusters in mass and redshift accord- 
ing to eq. © and eq. I|10(l . above the limiting mass dis- 
cussed above (see Fig.[(jJ). After excluding the nearby clus- 
ters at redshifts below 0.05, we obtain a catalog of ~ 1200 
clusters. 

We have then calculated the flux density ratios for var- 
ious frequency pairs and for different values of the param- 
eter a in the cosmological redshift-temperature law. We 
have used the bandwidths of the Planck instruments and 
calculated the observed SZ flux density assuming a top- 
hat response of the filters (R — 1 for v c — Av/2 < v < 
v c + /S.v/2 and R — elsewhere (see Table 1)). 

5.4. Error estimates 

Errors are simulated through use of the noise levels quoted 
in Table 1. As for the uncertainties on the flux density 



ratios, we have propagated the uncertainties on the indi- 
vidual flux densities. 

The random noise and the contribution from the cos- 
mic infrared backgrou nd and from unresolve d radiosources 
have been taken from lAghanim et alJ l|2004h and added in 
quadrature. Dust emission has been neglected as the levels 
of contamination are orders of magnitude lower than for 
the other contaminants mentioned above. 

Another contaminating source is the primary tempera- 
ture fluctuations of the CMB, which are important for ob- 
servations with a large beam such as that of Planck. The 
primordial anisotropics have the same spectral signature 
as the kinetic SZ effect. Promising techniques to separate 
the primary CMB fluctuations from the kinetic SZ effect 
have been developed, for example by taking advantage of 
the spatial correlation between the thermal and the kinetic 
effects (see lForni fc Aghaniml2004l and references therein) . 
In the present study, we have assumed that the contam- 
ination from the CMB could be removed. Observations 
with a smaller telescope beam close to the average size of 
the cluster (around 1 arcmin) wouldn't suffer significantly 
from contamination from the CMB. It has also been shown 
that such observations would pro yide an optimum surve y 
yield for the detection of clusters llBattve fc Wellerll2004l) . 

Although we assume precise and unbiased knowledge 
of the electron temperatures in the catalog, it should be 
pointed out that the results are not very sensitive to the 
cluster temperatures, as shown in Fig. 4. Adding an un- 
certainty of up to 20% on the cluster temperatures doesn't 
affect the final result significantly. Errors due to peculiar 
velocities have to be modelled more carefully, however; not 
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Fig. 8. Kolmogorov-Smirnov test: The significance level 
p is shown as a function of a. p is the probability that 
the two data sets are drawn from the same distribution. 
Stars show the results for the 100/143 GHz frequency pair, 
triangles for the 143/353 GHz pair. The 100/143 GHz pair 
appears as the better one, as the probability that the a = 
distribution be the same as a one with o / drops most 
rapidly with increasing a. Each data set contains 1000 
clusters. 



1. First assuming a — 0, we do a \ 2 analysis to fit the 
peculiar velocities to the simulated flux density ratios 
at two frequencies, namely 143 and 353 GHz, which is 
the frequency pair which is most sensitive to peculiar 
velocities. We do this for each cluster. 

2. Using those velocities, we then do a \ 2 analysis to fit 
the whole catalog to a, now using the frequency pair 
100/143 GHz, which is most sensitive to a. 

3. Taking the parameter a obtained at step 2, we iterate 
and fit to the peculiar velocities, and then fit to a new 
value of a. 

We iterate until the method converges. Figure [3] shows 
clearly that it does: Here we have deliberately used a large 
a: 0.4. After a few iterations, the method converges. We 
have also examined the distribution of the peculiar veloc- 
ities at each iteration. At the first iteration (when a is 
taken equal to zero), the distribution of peculiar veloci- 
ties is clearly non-Gaussian, as indicated by its skewness. 
With each integration, the distribution of velocities con- 
verges toward the original one; the parameter a also tends 
toward its original value. In order to test the robustness 
of the method, we repeated 200 times a simulation with 
a = 0. The distribution of the inferred a parameters fol- 
lows a Gaussian, with a dispersion around 10~ 2 . 



only because of the limited precision with which these can 
be extracted from the SZ data, but also because A/1/A/2 
is more sensitive to v r than to kT e . 

We use a separate frequency ratio to fit to v r for each 
individual cluster. Error bars on each peculiar velocity are 
converted into error bars in the ratio used for the CMB 
temperature by simply calculating this ratio at the end- 
points of the 68% confidence interval of the peculiar veloc- 
ities. These systematic errrors are then added in quadra- 
ture to the other errors. 

5.5. Analysis 

As a first step, we have compared simulated SZ obser- 
vations of catalogs with different a. Figure shows the 
flux density ratio for the frequency pair 100/143 GHz as a 
function of redshift in two cases: first, in a simulated cat- 
alog where all clusters have zero peculiar velocity (left), 
and then where peculiar velocities have been assigned ac- 
cording to eq. 1|17|) . Clearly, the introduction of peculiar 
velocities confuses the picture and makes it more difficult 
to distinguish the distributions with different a. In order 
to compare quantitatively the distributions with the in- 
clusion of peculiar velocities, we have run a Kolmogorov- 
Smirnov test (Fig. [HJ . We see clearly that the distributions 
are statistically different for a sufficiently large a. 

This first test is encouraging. In practice, however, one 
would like to be able to determine the value of the a pa- 
rameter from one single simulation, since we have only one 
observable universe at our disposal! To do this, we have 
devised the following iterative method: 



6. Conclusion 

We have examined the possibility to constrain the cosmo- 
logical temperature-redshift law from multifrcquency SZ 
observations. The advantage of using flux density ratios 
is to eliminate a possible dependence on cluster param- 
eters, in particular the optical depth and, to some ex- 
tent, the cluster temperature. After several tests, we have 
selected the 100/143 GHz frequency pair as particularly 
sensitive to our parametrization of the T — z law, and the 
143/353 GHz pair as sensitive to the cluster peculiar ve- 
locities. The iterative method that we have devised makes 
it possible to recover the parameter a of the T — z law 
from SZ flux density ratios of a large number of clusters. 
Our simulated catalog contains ~ 1200 clusters. We have 
assumed a ACDM cosmology, a full-sky survey and de- 
tection at more than three times the expected sensitivity 
of the Planck instrument taking into acccount contamina- 
tion by the Cosmic Infrared Background and radiosources. 
Note, however, that the method is independent of the cos- 
mological parameters and could be applied to any catalog 
of clusters for which redshifts and, to some extent, tem- 
peratures are known. 

Most cosmological models predict a linear variation of 
the CMB temperature with redshift. Sensitive SZ observa- 
tions appear as a potentially useful tool to test the stan- 
dard law. The discovery of an alternative law would have 
profound implications on the cosmological model, imply- 
ing creation of energy in a manner that would still main- 
tain the black-body shape of the CMB spectrum at red- 
shift zero. 
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Fig. 9. Convergence of the method. Here a simulation was a = 0.4 was performed. The method clearly converges 
toward that value of a, as can also be seen from several diagnositics: the mean value of the peculiar velocity distribution 
and the skewness converge toward zero as a tends toward its original value. 
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